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Abstract 

The one-loop QCD effective charge a e /f for quark-quark scattering is derived by di- 
agrammatic resummation of the one-loop amplitude using an arbitary covariant gauge. 
Except for the particular choice of gauge parameter £ = —3, ay* is found to increase 
with increasing physical scale, Q, as InQ or ln 2 Q. For £ = —3, a e J* decreases with 
increasing Q and satisfies a renormalisation group equation. Also, except for the case 
£ = 19/9, convergence radii of geometric series are found to impose upper limits on Q. 
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1 Introduction 



Quark-quark scattering in next-to-leading order QCD has been calculated by several dif- 
ferent groups [1-4]. Coqueraux and De Rafael [1], calculated the one- loop corrections to the 
invariant amplitude in the Feynman gauge using an on-shell renormalisation scheme [5]. In 
Refs. [2-4] complete expressions for the squared invariant amplitude were given in the dimen- 
sional regularisation scheme [6]. The calculations presented in the present paper generalise 
those of Refs. [1-4] in two ways: 

(i) An arbitary covariant gauge is considered. 

(ii) The one-loop Ultra- Violet (UV) divergent loop and vertex diagrams are resummed to 
all orders in a s . 

(ii) yields a scale-dependent 'effective charge' as a factor in the invariant amplitude. The 
resummation is done, not by solving a renormalisation group equation (RGE), but by an 
exact sum of the relevant diagrams to give the QCD analogue of the Dyson sum of QED. 

The results are very surprising. The effective charge is gauge dependent at O(a^) and 
beyond, and except for one specific choice of gauge, does not display 'asymptotic freedom' 
but instead increases as \nQ or In 2 Q at large scales Q. Only for the same special choice 
of gauge, where contributions from vertex diagrams vanish, does the effective charge satisfy 
a RGE of the type that is valid in QED [7-11]. For this special choice of gauge (called 
'loop gauge') the effective charge decreases with increasing scale, but only to a fixed limit 
Ql, determined by the convergence radius of the geometric sum of gluon and fermion loops 
in the dressed gluon propagator. The measured value of a s suggests that Ql — 300 GeV, 
already of phenomenological importance at the Fermilab pp collider. 

The plan of the paper is as follows. In the next Section the UV divergent (before renor- 
malisation) one-loop corrections to the quark-quark scattering amplitude in Feynman gauge 
are derived from similar corrections to the quark-quark-gluon vertex given in Ref. [12]. The 
generalisation to an arbitary covariant gauge is made using results reported in Ref. [13]. In 
Section 3, the one-loop corrections are diagrammatically resummed to yield the QCD ana- 
logue of the Dyson sum of QED. In Section 4 the self-similarity and renormalisation group 
properties of the effective charge derived in Section 3 are discussed. In the final Section the 
classical proofs of the asymptotic freedom property of QCD in the literature are critically 
examined in the light of the results obtained in the previous Sections. Also briefly discussed 
are : (i) 'renormalons' (ii) the generalisation to higher loop order vacuum polarisation and 
vertex corrections and (iii) pinch technique and related calculations of proper self energy 
and vertex functions, both in QCD and in the Standard Electroweak Model. 

The present paper is the fifth in a series on fundamental physics aspects of perturba- 
tive QED and QCD: Ref. [15] discusses on-shell renormalisation and optimised perturbation 
theory in QED and QCD; Ref. [16] considers the role of non-vanishing fermion or gluon 
masses as physical regulators for diagrams with combined infra-red (IR) and ultra-violet 
(UV) divergences; the relation of such IR/UV divergent diagrams to the Lee-Naunberg [17] 
and Kinoshita [18] Theorems of QED and QCD is examined in Ref. [19]; convergence condi- 
tions for resummed physical amplitudes, as discussed in the present paper for quark-quark 
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scattering in QCD, are analysed for the analogous QED case of scattering of unequal mass 
charged fermions in Ref. [20]. 



2 The quark-quark scattering amplitude to 0(a 2 s ) in an 
arbitary covariant gauge 

The process considered is the scattering of two equal mass quarks through an angle of 
90° in their CM system. The lowest order diagram is shown in Fig. la. The four-vectors of 
the incoming (outgoing) quarks are p\,P2 (p3,Pi)- In this configuration the exchanged gluon 
has a virtuality t = —s/2 where: 

t = ( Pl -p 3 ) 2 = u=( Pl -p 4 ) 2 , (2.1) 
s = (P1+P2) 2 = (P3+P4) 2 - (2.2) 

Denoting the invariant amplitude corresponding to the diagram in Fig. la by M.^ (quark 
spin and colour indices are suppressed), then the 0(a 2 ) amplitude may be written as: 

M m =M<®Y,Ai, (2.3) 

i 

where the one-loop virtual corrections Ai are given by the diagrams shown in Figs, lb - li. 
Three topologically distinct types of diagrams occur: 

• vertex corrections as in Figs, lb, lc and the two similar diagrams given by the exchange 
13 <-> 24; 

• loop corrections (Figs, ld-lg); 

• box diagrams (Figs, lh, li). 

The leading logarithmic corrections, Ai, corresponding to the contributions of diagrams that 
are UV divergent before renormalisation, have been derived in a straightforward fashion 
from the corrections to the quark-quark-gluon vertex presented in Ref. [12]. The results, in 
Feynman gauge, derived using dimensional regularisation [6] are presented in Table l 1 . Here 

In two previous versions of the present paper [21] the one- loop corrections to the quark-quark scattering 
amplitude were taken from Table 1 of Ref. [1] . This calculation was performed using on-shell normalisation [5] 
in which explicit quark and gluon mass parameters were introduced. Comparison of the results of Ref. [1] 
with those of Refs. [12] and [13] and those for the related QED Bhabha scattering process [22] have revealed 
a number of inconsistencies with the results of Ref. [1] . (i) the ' Vertex(a + b) ' contribution shown is simply 
the QED result modified by a single multiplicative colour factor. The QCD-specific UV divergent term 
found in Refs. [12] and [13] and presented in the first row of Table 1 in the present paper is absent, (ii) 
The 'Three Gluon(i+j)' contribution should, according to Ref. [12], be 3y/2, not —y. (iii) Adding up the 
singly logarithmic terms of Table 1 of Ref. [1] yields, as the coefficient of the logarithm of the physical scale, 
/3o, but the contribution of 'Vertex(a and b)' actually corresponds, not to a UV divergent term, but an IR 
divergent one, that is cancelled on adding the contributions of diagrams with real gluon radiation from the 
external quark lines. Thus Eqs. (2.4) and (2.5) of the previous versions of the present paper were incorrect 
and the appearence of f3o in Eq. (2.12) was fortuitous. However, Eqs. (2.5) and (2.6) of the present revised 
and corrected version are identical to Eqs. (2.15) and (2.16) of the previous versions, so that none of the 
subsequent discussion or conclusions are affected by the corrections. 
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Q = y/^i = ys/2 is the physical scale, and /i is the renormalistion subtraction scale. Ca is 
the usual QCD colour factor (Ca = number of colours = 3) while n/ is the number of light 
quark flavours contributing to the vacuum polarisation loops in Fig. Id. is the square of 
the renormalised on-shell strong coupling constant at the scale \x. The corrections Ai—A^ 
may be combined to yield a term proportional to the first coefficient in the perturbation 
series in a s of the beta function of QCD. Denoting the quark loop + gluon loop + ghost 
loop correction, A3 + A4, by and the vertex correction, A\ + A 2 , by V M , and generalising 
to an arbitary covariant gauge specified by the parameter £, in which the gluon propagator 
is written as : 



Q 



<r-(i-o 



(2.4) 



yields the explicit gauge dependence 2 of V M and C,^ [13] 



v(0 



71 4 V 



Q. 



7T 



■-(--0 + ^ 



1 fi\ 



(2.5) 
(2.6) 



Adding Eqs. (2.5) and (2.6): 



< 11 n/ 



71 Z 6 fl TX fl 



where (3 is the first coefficient of the QCD (3 function [ a s = a s /ir] : 

-[3 a 2 s + ... 



da s 



To 0(a 2 s ), af may be identified with the solution, af GE (Q), of the RGE (2.8): 
m lo = mW+M™ =M m (l + V i (£) + 



O 

<[l--^/3 l„(^)] 

7T /i 



.M^af^Q) 



Ots 



+ o (K) 2 ) , 



where 



a 



RGE 



(Q) = «? = 



1 + f /3 ln©' 



(2.7) 



(2.8) 



(2.9) 



(2.10) 



.M LO denotes the amplitude including Leading Order vertex and vacuum polarisation correc- 
tions. In the following section the amplitude Ai^, in which these corrections are summed 
to all orders in a s , is derived. 

It can be seen that, at one loop, the gauge invariant result (2.7) is obtained on adding 
the UV divergent contributions of diagrams (b)-(g) in Fig. 1. The box diagrams (h) and (i) 



2 This gauge dependence of one- loop vertex and vacuum polarisation insertions containing triple gauge 
boson couplings is universal for all non-abelian gauge theories. See, for example, Ref. [14] Eqs. 12.114 and 
12.122. 
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and those obtained by exchanging the internal quark and gluon propagators, form, together 
with the diagrams where a single gluon is radiated from one of the external quark lines of 
Fig. la, another gauge invariant set. Indeed the contribution to the cross section of this set of 
diagrams, which involve only abelian quark-gluon couplings, is obtained by multiplying the 
result for the analogous QED t-channel Bhabha scattering process [22] by the appropriate 
QCD colour factor: C 2 (N C ) = (iV c 2 - 1)/2N C = 4/3. The contribution of the box diagrams 
alone is UV finite, but IR divergent for the case of massless gluons. This IR divergence is 
cancelled on adding the contribution of the real gluon radiation diagrams. The contributions 
of the box diagrams and the real radiation diagrams are however, separately, gauge invariant 3 . 
Since the box diagrams are UV finite, and so do not require renormalisation, they do not 
contribute to the QCD effective charge. All the terms that contribute to the latter at one 
loop are presented in Table 1; they are UV divergent before renormalisation. 



3 The resummed quark-quark scattering amplitude 



The topographical structures 4 of the diagrams that modify the gluon propagator in the 
quark-quark scattering amplitude at 0(a 2 ), 0(a 3 ) and 0(a 4 ) are shown in Figs. 2a,b,c 
respectively. V and L denote vertex and loop (vacuum polarisation) contributions: 

V = Ui + U 2 , (3.1) 

n f 

L = J2^ + Gi + G 2 + G 3 . (3.2) 
i=i 

V u V 2 correspond to the diagrams in Fig. lb, lc ; £?4 F t to Fig. Id and G u G 2 , G 3 to 
Figs. le,f,g. In the case that only one vertex insertion occurs there is a factor 2 for the two 
ends of the gluon propagator. Since the propagator has only two ends the vertex corrections 
are never higher than quadratic in the perturbation series for the amplitude. Although the 
topographical structure of diagrams containing vertex corrections is different at 0(a^) than 
at 0(a 2 ) it remains the same at all higher orders. The all-orders resummed amplitude is: 

M (oo) = M {0) [1 + 2V + L 

+ V 2 + 2VL + L 2 
+V 2 L + 2VL 2 + L 3 + ... 

= M {0) [l + L + L 2 + ... + V 2 (1 + L + L 2 + .. 

+2V(1 + L + L 2 + ... 

= M^_+V)\ ,3.3, 



3 See Ref. [23] for a complete diagrammatic discussion of gauge cancellations in the one-loop corrected 
quark-quark scattering amplitude. 

4 The use of the word 'topographical structure' to indicate a particular disposition of vertex and (possibly 
resummed) self-energy insertions in a diagram is deliberate. Diagrams with internal lines in the vertex or 
self-energy insertions have a different topology but the same topography as the one-loop diagrams shown in 
Fig. 2. 
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In an arbitary covariant gauge with momentum subtraction at scale /i, and in leading loga- 
rithmic approximation, V = V M (£)/2, L = so that 



_A/f(oo) = 



as 



1 - />(£) 



(3.4) 



leading to the following expression, in leading logarithmic approximation, for the resummed 
one-loop effective charge: 



al"(Q) = < 





l-lf(3 + 0ln(f); 


2 


1 + Si 


13 - 3£ - 


ln(2) 



(3.5) 



The conventional one-loop QCD running coupling constant Eq. (2.10) is recovered only for 
the special choice of gauge parameter £ = = —3 ('loop gauge'). Only in this case does 
a e Jf {Q) decrease monotonically with increasing Q. For any other choice of gauge a e /f{Q) 
does not show 'asymptotic freedom' as Q — > oo but instead increases as InQ when £ 7^ £y 
and as In 2 Q when £ = £y. The choice £ = £y ('vertex gauge') where 

£ y = (39-4n/)/9 

corresponds to a vanishing coefficient of ln(Q/yu) in the denominator of Eq. (3.5). As dis- 
cussed in more detail below, only in loop gauge is the equation for the effective charge 'self 
similar' like the effective charge in QED or the solution (2.10) of the RGE Eq. (2.8). Even 
in loop gauge the effective charge of Eq. (3.5) does not decrease without limit as Q — > 00. 
The maximum possible scale Ql (Landau scale) is determined by the convergence proper- 
ties of the geometric sum that yields the denominator of Eq. (3.5). The geometric series is 
convergent provided that \L\ < 1 [24]. This implies that: 



a 



4:71 



3£ 



The corresponding Landau scale is then: 



Qi 



/xexp 



47T 



< 1. 



(3.6) 



a£[13-3£ 



4nq 
3 J 



and 



For £ = -3: 



a 



eff 



a' 



(Ql) = -rf 



3< 
~8r7 



(3 



Oln(%) 



< ff (Q L ) = f , 



(3.7) 



(3.8) 



(3.9) 



so the convergence condition (3.6) implies that in this case a e J* cannot evolve down by more 
than a factor of | of its initial value before Eq. (3.5) diverges — there is no 'asymptotic 
freedom'. The convergence conditions for a geometric series with a negative common ratio, 
such as that which generates the QCD running coupling constant, are derived below in an 
Appendix. Unlike in the case of the Landau pole [25] of QED, with a positive common ratio, 
r = \r\, where the running coupling constant is proportional to 1/(1 — |r|), it is not obvious, 
by inspection, that the QCD perturbation series is not equal to 1/(1 + |r|) when |r| > 1. 
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However (see the Appendix) this is quite clear from the exact formula, Eq. (Al), giving the 
sum for any finite number of terms. Physicists have the right to be as free as possible in 
making conjectures in their attempts to describe nature in the simplest way possible, but 
not to make ones that are contrary to mathematical laws. 

Numerical values of Q L for /i = 5 GeV, = 0.2 and rif — 5 (corresponding, approxi- 
mately, to the experimental value of a£) and four different choices of gauge parameter are 
presented in Table 2. For loop gauge (£ = —3) Ql is only 300 GeV, with phenomenological 
consequences perhaps already at the Fermilab pp collider, but certainly at the future LHC 
pp collider. It may be noted that, in all cases except vertex gauge, Ql lies well below the 
Grand Unification (GUT) scale of ~ 10 15 GeV. Since for any gauge choice except £ = —3, 
oi e Jf diverges as InQ or In 2 Q at large Q, there can be no 'unification' [26, 27] of the strong 
and electromagnetic interactions at large scales Q, for any choice of gauge parameter, at 
least if the running strong coupling constant is identified with an effective charge such as 
that in Eq. (3.5). In fact all studies, to date, of Grand Unification have implicitly used loop 
gauge where the maximum scale Ql is only ~ 300 GeV. In vertex gauge, since all vacuum 
polarisation contributions vanish, there is no convergence limitation on Q in Eq. (3.5). 

The behaviour of a e J* at large values of Q depends on the value of the Landau scale and 
value, Qo, of Q at which the numerator of Eq. (3.5) vanishes: 



Qo = H exp 



8tt 



_3a£(3 + 0. ' 

The value of £, £ 0? at which Q L and Q are equal is given by (3.7) and (3.10) as: 



(3.10) 



eo(n/) = ifl7-^) (3.11) 



9 V 3 

For rif = 5, £o = 0.4074.... When £ < £ , (e.g. Loop and Landau gauges) and Ql < Qo, then 
a e /- f diverges before vanishing and so decreases monotonically with increasing Q within its 
convergence radius. For £ > £ (e.g. Feynman and Vertex gauges) Ql > Qo- In this case 
a e Jf first decreases with increasing Q, reaching a minimum value at zero, and then increases 
with increasing Q up to Ql- In vertex gauge Ql is infinite. Numerical values of Q for the 
same parameter choices as in Table 2 are presented in Table 3. The scale Qo is defined (i.e. 
Ql > Qo-) only for Feynman gauge (Q = 176 TeV) and Vertex gauge {Qo = 18.1 TeV). 

The diagrams shown in Fig. 2 each contain only a single 'dressed' virtual gluon. It is 
also possible to consider cases in which the subsitution of the series of diagrams illustrated 
in Fig. 2, and summed in Eq. (3.3), is made in a diagram containing more than one virtual 
gluon (for example the box diagrams shown in Fig. lh and li). This will result in a further 
higher order correction to the quark-quark scattering process. It is clear, however, that 
the manifest gauge dependence of the resummed one-loop contribution will be unaffected 
by the presence of other virtual gluons ('dressed' or not) in the diagram. Every virtual 
gluon has only two ends, and so the resummed vertex correction is at most quadratic. Then 
no cancellation is possible of the gauge dependent pieces of multiple vacuum polarisation 
insertions. However, as further discussed in Section 5 below, the box diagrams themselves 
can be resummed. This typically results in a Sudakov-like double logarithm, not a geometric 
series as found for resummed vacuum polarisation diagrams. 
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4 Self-Similarity Properties of a e Jf and the Renormal- 
isation Group 



Introducing the abbreviated notation: 

3 



«(0 = g(3 + 

a e /f(Q)/7r = a Q 
Eq. (3.5) may be written as: 



• '(f) = i 
i = ^ , 

q M [l - a^(QA] 2 
l + a M Z(0A 

With £ = -3, u(-3) = 0, Eq. (4.1) becomes 



13-3^ 



4ra, 



A = ln(Q//z) 



1 + a M /( 



-3) A' 



(4.1) 



(4.2) 



Eq. (4.2) is the solution of a one-loop RGE similar to that for the QED effective charge [10, 
11]: 

Q da Q U o \ a 

-l{-3)a Q = -fJ a Q , 



a Q dQ 



where 



For a gauge choice such that v ^ the partial differential equation satisfied by aq is: 



Q da Q 
a Q dQ 



1(0 



o M v(0A) : 



+ 2 



(1 - a^(OA) 3 



(4.3) 



(4.4) 



(4.5) 



so that, in this case, the effective charge oq does not satisfy the RGE (4.3). Expanding in 
powers of a M on the right hand side of Eq. (4.5) gives: 



Q da Q 
a Q dQ 



= -a Q [l(t;) + 2v(0 + O(a fl )} 
= ~a Q [A, + 0(a M )] . 



(4.6) 
(4.7) 



Thus, neglecting terms of 0(aQa M ) ~ O(a^) so that only the first term in the QCD per- 
turbation series is retained, oq — > aq and it can be seen that Oq , for an arbitary gauge 
choice, satisfies the same partial differential equation as the one-loop all orders resummed 
ag in loop gauge. The relation of this result to previous derivations of the QCD running 
coupling constant, where it has generally been conjectured that a gauge invariant result is 
obtained to all orders in perturbation theory, is discussed in the following Section. 

The equation (4.2) is self-similar in the sense that , for any values of fi and Q the equation 
defined by the exchange /i <-> Q is identical to the original equation. A consequence of this 
symmetry property is that Oq in Eq. (4.2) is independent of fi (with the important caveat 
that, since the denominator is the sum to infinity of a geometric series, \x must be such that 
\a^l(— 3)A| < 1), and that in the equation with /i <-> Q, is independent of Q. This is the 
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mathematical basis of the Renormalisation Group [7, 8, 9]. Such a universal self-similarity 
property is not, however, shared by Eq. (4.1) when v(£) ^ 0. 

Eq. (4.1) is self-similar under the exchange li <-> Q provided that the equation: 

a Q [l + a Q v(QX] 2 

and Eq. (4.1) are both valid. Simultaneous solution of Eqs. (4.1), (4.8) in the case that 
v (0 leads to a quadratic equation for a M with the solution: 



Real solutions of Eq. (4.9) exist provided that either 



>a Q > — - (4.10) 



(both factors under the square root positive) or 

M^ > ± > _!^ (4 . n) 

2 a Q 2 v ; 

(both factors under the square root negative). 

The solution for a M , Eq. (4.9), is independent of /(£). Thus for fixed values of oq, Q, li, 
v(£) 7^ 0, satisfying the conditions (4.10) or (4.11) there are, in general, two values of such 
that Eq. (4.1) is self-similar. This value of a M has however no relation to the effective charge 
at the scale li given by Eq. (4.1) when Q — li. Choosing Q = 90 GeV and Landau gauge then 
(see Fig. 3) a Q = 0.095/tt = 0.0302. The choice li = 5 GeV gives 2/(u(0)A) = 1.933. The a Q 
terms under the square root of Eq. (4.9) may then be neglected, leading to the solutions: 

a + — + — = 1 9648 

a » ~ v(Q)X + 2 i " yb48 

a~ ~ = 0.0152 
2 

to be compared with the physical value: 

= a eff (5GeV)/n = 0.2/tt = 0.064 



Unlike for the special case v — 0, the first derivative of aQ in general depends on the 
scale li. For Q = ii the derivative is a negative constant fixed by the first coefficient in the 
perturbation series for the beta function as in Eq. (4.3). For any other choice of li when 
v 7^ the derivative varies with Q and li according to Eq. (4.5). Thus, if the effective charge 
is parametrised in terms of a M , cy at the reference scale Qr\ 

Qr ~ " [1 + aM)M ~ ' [1 + V'(0Ay ' ( ] 
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X R EE HQR/V) > = ^(Qfl/A*') 

then for Q 7^ Q_r and i>(£) 7^ the effective charge clq predicted by the formula containing 
fx (the second member of Eq. (4.12)) will be different to that predicted by that containing 
fi f (the third member of Eq. (4.12)), so that the value of czq depends upon the choice of 
renormalisation scale — it is no longer invariant. Renormalisation group invariance with 
respect to the choice of the scale \x is therefore not respected unless v(£) = 0, i.e. £ = —3. 

5 Discussion 

In the original derivations of the 'asymptotic freedom' property of QCD [28, 29] no calcu- 
lations were performed beyond the lowest non-trivial order, O(o%), and no actual amplitudes 
for physical processes were considered. It was conjectured (without any check by direct dia- 
grammatic calculation) that the QCD running coupling constant (RCC) could, in general (for 
any choice of gauge) be identified with the solution of the differential equation (2.8). The one- 
loop QCD beta function was calculated by considering the Callan-Symanzik [10, 11] equation 
for an irreducible n-point function (typically the gluon-quark-quark vertex or the triple gluon 
vertex). Calculation of the anomalous dimensions of the quark and gluon fields then yields 
the (gauge invariant) expression for the first beta function coefficient (3q in Eq. (2.8) above. 
An analogous result is obtained above by considering the unresummed one-loop correction 
to the physical quark-quark scattering amplitude. The true high order behaviour, however, 
corresponds to the sum of all the possible amplitudes for the process of interest. For this the 
actual topographical structure of the diagrams contributing to the amplitude must be prop- 
erly taken into account. Renormalisation scale invariance of the RCC, and the asymptotic 
freedom property, for an arbitary choice of covariant gauge, are not confirmed by the dia- 
grammatic calculation of higher order corrections in the case of the quark-quark scattering 
amplitude considered in this paper. 

Gauge dependence of the RCC in QCD has been considered previously in the liter- 
ature, but usually as an effect only at the two-loop level and at higher orders. It was 
pointed out that, in the case when the bare parameters of the theory are held fixed, the 
gauge parameter becomes scale dependent, and for certain momentum subtraction renor- 
malisation schemes, the second coefficient of the beta function is both scheme and gauge 
dependent [30]. For an arbitary covariant gauge specified by the fixed parameter £, as in 
the one-loop discussion above, the second beta function coefficient is, however, both gauge 
and renormalisation scheme invariant. Assuming that the RCC, mass and gauge parameter 
each satisfy renormalisation group equations similar to (2.8), and solving the coupled system 
of differential equations, solutions were found for the RCC that strongly depended on the 
initial conditions imposed on the running gauge parameter [30]. These solutions exhibit ei- 
ther asymptotic freedom-like behaviour or increase with increasing scale until an ultraviolet 
fixed point is reached [31, 32]. As commonly done in the literature, the RCC was treated as 
an independent mathematical object, without reference to any actual physical process, and 
the renormalisation group equations were assumed to hold without specific diagrammatic 
justification. It is shown above that, if the RCC is identified with the effective charge of the 
quark-quark scattering amplitude, the one-loop renormalisation group equation 5 holds only 

5 In Rcfs. [31, 32] the one-loop RGE was assumed to be gauge independent and given by the conventional 
formula (2.10). 
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for the specific gauge choice £ = — 3. The gauge parameter can then neither vary nor satisfy 
a RGE. 



When quark mass effects are taken into account, the one- loop beta function coefficient 
is also gauge dependent, and has a value which depends on the particular n-point Green's 
function considered for its derivation. The mass dependent corrections to the triple gluon 
vertex [33] and the gluon-ghost-ghost vertex [34] are different. A detailed discussion may be 
found in Ref. [35]. A corollary is that the RCC in physical amplitudes is both gauge and 
process dependent at physical scales where quark mass effects (other than those contained 
in the asymptotically dominant logarithmic terms) are important. 

The effective charge (3.7) has been calculated here for the simple case of quark-quark 
scattering with a unique physical scale Q = t. In this case the direct physical interpre- 
tation as the strength of the interaction between two currents varying as a function of their 
separation (~ 1/Q) is particularly transparent. However, since every dressed propagator 
has just two ends, similar expressions for the RCC (in general a function of some running 
loop 4-momentum k) are expected in all physical amplitudes containing virtual gluon lines. 
Two examples are shown in Fig. 4. Fig. 4a shows the three topographically distinct classes 
of diagrams that contribute to the anomalous magnetic moment of a heavy quark at 0(a^). 
In Fig. 4b the same classes of diagrams are shown for the four-loop photon proper self 
energy function due to radiatively corrected quark vacuum polarisation loops. As for the 
quark-quark scattering case the same topographical structures (giving at most a quadratic 
dependence on the vertex corrections) is found at all higher orders in the 'dressed gluon 
propagator'. The diagrams of Fig. 4b are related via the optical theorem and analytical 
continuation to the process: 

e + e~ — > 7* — > qq + X 

where X denotes g, gg or q'q'. 

There has been considerable recent interest in the structure, in high orders of pertur- 
bation theory, of diagrams containing chains of vacuum polarisation loops in internal gluon 
propagators (for example the generalisation to higher orders of the 0(af) diagrams with 
two vacuum polarisation loops shown in Fig. 4). The anstatz used for these so-called 'renor- 
malon chains' [36] is to replace ri//3 in a calculation considering only rif different flavours 
of fermion vacuum polarisation loops with nf/3 — 11/2 = —(3q. This is clearly a good ap- 
proximation in the limit nj — > oo. As inspection of Fig. 4 shows, however, this will result, 
in any gauge in which the vertex correction is non-vanishing, in a miscounting of the con- 
tribution of the latter, which are included at order n in terms of the form (3q, but actually 
should never appear at higher order than quadratic in the perturbation series. With how- 
ever the gauge choice £ = —3 (loop gauge) all vertex corrections vanish and the renormalon 
chains are correctly given by the above replacement. This gauge choice is, in any case, the 
one universally (although tacitly) made in all phenomenological applications of the QCD 
running coupling constant, where renormalisation group invariance is assumed. Thus the 
'Naive Non-Abelianization' (NNA) ansatz [37] or 'Large /3q Limit' that is typically assumed 
in phenomenological studies of renormalon effects is a correct one only in loop gauge. The 
remark that the choice of gauge £ = —3 implies the vanishing of all vertex corrections, has 
previously been made in the context of two-loop corrections to heavy quark production in 
e + e~ annihilation [41]. 
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An instructive example of the inconsistent treatment of Feynman diagrams by the NNA 
anstatz as well as a illustration of the gauge dependence of fixed order perturbative QCD cal- 
culations beyond next-to-leading-order (NLO) is provided by the analysis of the moments of 
non-singlet anomalous dimensions of deep-inelastic nucleon structure functions in Ref. [38]. 
In this work the effect of insertion of an arbitary number of quark or gluon vacuum polari- 
sation insertions in the virtual gluon propagator of a forward Compton scattering amplitude 
(or the equivalent diagrams in the Operator Product Expansion, as in Fig. 1 of Ref. [38]) is 
considered. An example of an 0(a s ) diagram that contributes to this amplitude is shown 
in Fig. 5a. One and two vacuum polarisation loop insertions as considered in Ref. [38] are 
shown in the fourth topographical diagrams of Fig. 5b and 5d respectively. The predic- 
tions of this procedure for the loop (£ = —3) and Landau (£ = 0) gauges were compared 
to the exact massless next-to-next-to-leading-order (NNLO) calculation in Feynman gauge 6 
(£ = 1) of Ref. [42]. Table 4, extracted from Table 1 of Ref. [38], shows the results of this 
comparison for the n = 2 moment at NLO, as in Fig. 5b, and NNLO, as in Fig. 5c and 5d. 
In the calculations of Ref. [42] the contributions of all 353 Feynman diagrams contributing 
to the anomalous dimensions up to NNLO were evaluated, whereas in Ref. [38] only quark 
and gluon loop vacuum polarisation insertions in Fig. 5a and the other Compton scattering 
diagrams with one virtual gluon line (see for example Fig. 1 of Ref. [43]) were considered. 
Even so, agreement is found with the Loop gauge renormalon calculation at the 20% level 
at NLO for the n = 2 moment. Even better agreement is found for higher moments — it is 
good at 4.7 % for the n = 6 moment. However, at NNLO no agreement is found; indeed 
predictions of the 'renormalon dominated' approximation with either choice of gauge have 
even a different sign to the exact NNLO calculation. The reasonable agreement at NLO 
between the two calculations when loop gauge is employed in Ref. [38] can be understood by 
inspection of Fig. 5b. This choice of gauge is equivalent, at NLO, to performing a calculation 
in an arbitary covariant gauge, in which the non-abelian one-loop corrections to the quark- 
gluon coupling are included as well as vacuum polarisation insertions. All the NLO vertex 
and vacuum polarisation corrections are included in the topographical diagrams of Fig. 5b. 
The sum of these contributions is gauge invariant. The situation is quite different at NNLO. 
The subset of NNLO diagrams shown in Fig. 5c, obtained from those of Fig. 5b by inserting 
an additional virtual gluon between the incoming and outgoing quark lines also gives a gauge 
invariant result and is correctly described in loop gauge. This is not the case for the diagrams 
of Fig. 5d — the corresponding sum of the diagrams is of the form V 2 + 2VL + L 2 , which 
is manifestly gauge dependent. Other NNLO contributions with the same colour factors as 
in Table 4 arise from irreducible two-loop vertex and vacuum polarisation diagrams. The 
topographical pattern is the same as in Fig. 5b with the one-loop insertions replaced by ir- 
reducible two-loop insertions. As discussed in more detail below, if these two-loop insertions 
satisfy a Ward identity similar to that respected by the one-loop insertions, this contribution 
will also be gauge invariant. No possibility exists however to remove the gauge dependence 
of the contribution of the diagrams shown in Fig. 5d. This explains the breakdown of gauge 
invariance shown by the CfC\ NNLO entries of Table 4. 

6 The formulae given in Ref. [42] are actually in Landau gauge. However, it is stated in the paper that 
for low order moments n = 2, 4 "...the diagrams were run with a gauge parameter £ in the gluon progagator 
Qiiv — £,qfj,qu / q 2 ■" The assumed value (or values) of their parameter £ (1 — £ in the notation of the present 
paper) were not stated and for the moments n — 6, 8 the gauge parameter was not included (presumably it 
was set to zero in the calculations) which corresponds to choosing Feynman gauge. It is assumed here that 
results compatible with this choice of gauge were also obtained for the n — 2, 4 moments. 
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In an attempt to give a diagrammatic justification of the NNA ansatz, Beneke in Ref. [36] 
considered NLO loop corrections as in Fig. 2a or Fig. 5b above, for the case of quark pair 
production by a vector current. Performing the calculation in Landau gauge it was found 
that the beta function of the corresponding QED calculation was replaced by the one-loop 
QCD beta function. If this calculation had been performed in an arbitary covariant gauge, 
the gauge invariance of the QCD beta function would have been demonstrated. This in 
no way justifies the general use of the NNA ansatz since manifest gauge dependence as in 
Fig. 2b or Fig. 5d first occurs at NNLO where the gauge dependence of the V 2 term is not 
cancelled by NNLO vacuum polarisation contributions. 

In connection with the work presented in the present paper it is important to notice 
that the renormalon singularities arise due to loop integrals over the virtuality of an internal 
photon or gluon line in a Feynman diagram, as in Figs. 4 and 5. Although the renormalon 
singularities are related to singular IR or UV behaviour of the running coupling constant — as 
discussed by Lautrup in Ref. [36] , in QED it is the UV Landau singularity — the renormalon 
singularity occurs for arbitary values of the external physical scale due to the infinite range 
of the internal loop momentum 7 . In contrast, in the simple case of quark-quark scattering 
discussed in the present paper, there is no integration over the virtuality of the gluon line 
in which the vertex and loop corrections are inserted and the scale in the running coupling 
constant is identical to the physical scale of the problem. The diagramatic analysis shown 
in Fig. 2 is therefore much simpler and the breakdown of gauge invariance appears as soon 
as the NNLO contributions of Fig. 2b are evaluated. 

The classification of diagrams as in Figs. 2, 4 and 5, according to the categories (in an 
obvious notation) L n , VL n , V 2 L n , will remain when L, V are calculated with an arbitary 
number of internal lines. The global structure of Eq. (3.5) will then remain the same for 
calculations including an arbitary number of loops, though additional non-leading terms in 
In Q will result from integration over internal loops in the basic one-loop vacuum polarisation 
and vertex diagrams shown in Fig. 1. 

As in Ref. [1], only the leading logarithmic terms in the one-loop correction (and hence 
in the resummed effective charge) have been taken into account in the above discussion. 
Constant terms in V and L have been neglected. For a general renormalisation scheme 
however, constant gauge and renormalisation scheme dependent terms also occur, so that 
Eq. (4.1) is replaced by the expression: 

a °- a " 1 + «,[«{)A + „({)] ■ (!U) 



In the MS scheme [23]: 



/>x lQrif 97 3, 3 ~ 

c,n = — *- e — £ 2 

vs; 9 12 8 16 



7 In a discussion of renormalons in a review talk by S. Forte [44] the following important statement can be 
found: 'If the series had alternating signs the singularity (ultraviolet renormalon) would not be on the path 
of the integral' (in the Borel transform) 'but the integral would still run outside the radius of convergence 
of the series; we will not discuss this any further.' This is the only place in the literature, to the present 
writer's knowledge, where the limited domain of convergence, in the UV limit, of the QCD RCC (discussed 
in detail in the present paper) is mentioned. The perturbation series corresponding to the RCC of QCD 
indeed has 'alternating signs'. 



12 



By a suitable scale choice \i = jj! and with £ = 



—3 Eq. (5.1) may be written as: 



aQ = l + v [/(-3)lnf + q(-3)]' (5-2) 

So in this (non asymptotic) case even in loop gauge the effective charge does not corre- 
spond exactly to the solution (2.10) of the one-loop RGE (2.8). Numerically l(— 3) = 3.833, 
q(— 3) = —3.090 for ny = 5. Thus, in the MS scheme in loop gauge, the resummed one- 
loop invariant charge is only asymptotically a renormalisation group invariant when constant 
terms in the one-loop correction are retained. 

Following the observation that the UV divergent parts of the vertex corrections in Figs, lb 
and lc may be associated with related diagrams in which the virtual quark propagators 
are shrunk to a point (or 'pinched') it was suggested [45, 46] to redefine a gluon proper 
self energy function by adding to the contributions of Figs, ld-lg, that of the pinched 
vertex diagrams. At one-loop order the resulting gluon proper self energy function is then 
gauge invariant. It was then (incorrectly) stated that a gauge invariant resummed gluon 
propagator may be trivially derived from the one-loop result (for example, Eq. (2.19) of 
Ref. [46]). In fact it is easy to show, quite generally, that if the one-loop corrected quark- 
quark scattering amplitude is gauge invariant (the correct initial assumption of the pinch 
technique calculations of Refs. [45, 46]) then resummed amplitudes at all higher orders must 
be gauge dependent. Introducing the gauge invariant one- loop quantity: 

B = L + 2V. (5.3) 

The resummed amplitude at 0(a™ +2 ) may then be written for n > 1 (see Fig. 2 and Eq. (3.3)) 
as: 

M (n+2) = M {o)( L + V ) 2 L n -\ (5.4) 

Expressing J^A^ n+2 ^> in terms of the gauge invariant quantity B and the gauge dependent 
quantity L(£) gives: 

M (n+2) = M (0)}_( B + L(^)2 jL (^n-l > (5,5) 

which is manifestly gauge dependent. The Dyson sum in Eq. (2.19) of Ref. [34] correctly 
describes the all orders resummed amplitude, not for an arbitary gauge parameter £, but 
only for the special choice £ = — 3 when V(£) = 0, B = L(—3) and 

M (n+2) = M (0) B n+l ^ = _ 3 y ^ 

Clearly, the above argument for manifest gauge dependence, shown to be valid at the re- 
summed one-loop level must also hold at arbitary loop order if the vertex and self-energy 
insertions satisfy a generalised Ward identity giving, at each order of perturbation theory, 
a condition such as (5.3). It has been shown [47], by the application of background field 
techniques, that Ward identities relating vertex and self-energy contributions may indeed 
be derived that are valid to all orders in perturbation theory. The gauge independence of 
the Ward identity at each fixed order then necessarily implies gauge dependence when the 
corresponding vertex and self-energy diagrams are resummed. 

The manifest gauge dependence of the quark-quark scattering amplitude found, by direct 
calculation, in this paper, is, apparently, in contradiction with formal proofs [48, 49] of the 
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gauge invariance of S-matrix elements in non-abelian gauge theories. It seems however, 
that what is actually proved in these papers is the gauge invariance, at all orders in pertur- 
bation theory, of generalised Ward identities. The consequences of resumming diagrams of 
fixed loop order, which as shown above, necessarily generates gauge dependence, were not 
considered. For example, in the standard reference [49], the Lagrangian from which all the 
Feynman rules of the theory is derived is introduced, and the change in this Lagrangian 
due to a change in the gauge parameter written down. It is then stated that the theory 
is gauge invariant if such a variation of the gauge parameter leaves S-matrix elements in- 
variant. There immediately follows the statement: 'We can formulate this condition' (i.e 
that the S'-matrix elements are invariant) 'in terms of a Ward identity that we have written 
in terms of the diagrams in Fig. 2'. The unproved and unjustified assertion is thus made 
that the gauge invariance of a Ward identity is equivalent to gauge invariance of S'-matrix 
elements. This will only be true of unresummed amplitudes at each loop order, not of the 
resummed amplitudes that, according to quantum mechanical superposition, must exist and 
are, indeed, essential to generate the RCC. In fact Ref. [49] establishes only the gauge invari- 
ance of Ward identities, and nothing else. As shown above, it is just the gauge invariance 
of the unresummed amplitudes that ensure the manifest gauge dependence of the resummed 
ones. Indeed an S-matrix element (even a formal, generic, one) appears nowhere among 
the equations of Ref. [49]. In Ref. [48] such a formal S-matrix element does appear, but its 
gauge invariance properties are derived directly from a Ward identity. No actual physical 
process, and no effect of resummation, is considered. 

By consideration of a sub-set of n-loop diagrams for the off-shell gluon-gluon scattering 
amplitude in a non-covariant gauge, it has been claimed [23] to demonstrate that the RCC 
of QCD is both gauge invariant and process independent, and that it may be identified 
with the solution (2.10) of the RGE (2.8). The n-loop diagrams considered are those that 
may be constructed as a formal 'product' of n + 1 tree level four-point functions. The 
diagrams contain both resummed one-loop gluon vacuum polarisation and vertex diagrams 
and a sub set of irreducible n-loop diagrams. This set of diagrams cannot, as claimed, be 
identified with the one-loop RCC in Eq. (2.10), which results solely from the resummation 
of one-loop (one particle irreducible) diagrams. At any order in the perturbation series 
these resummed one- loop diagrams contribute the leading powers of both In Q and rif. They 
give, in fact, the 'renormalon' contribution [36] (see above) that dominates the high order 
behaviour of the perturbation series. The irreducible n-loop (n > 1) diagrams of the sub-set 
considered in Ref [23] will contribute constant terms or non-leading powers of InQ, and 
therefore cannot be identified with terms in the diagrammatic expansion of the RCC in 
(2.10). Similarly, it has been conjectured (without explicit calculation) in Ref. [50] that the 
'missing' vertex contributions needed to make, say, Ai^ in Eq. (5.5) above, gauge invariant 
may be derived from 'pinch parts' of two-loop irreducible diagrams. This is not possible since 
the required 'missing' contributions contain the factor (a s \nQ) 2 , (see Eq. (3.5), whereas, as 
is well known, in both QED [51] and QCD [52] irreducible two-loop vacuum polarisation 
and vertex diagrams have, at most, next-to-leading logarithmic behaviour a 2 InQ. This 
is easily demonstrated by considering the two-loop solution of the renormalisation group 
equation for the effective charge. In QED, or in QCD in the £ = —3 gauge, Eq. (4.2) 
generalises to: 

°° = l + a M A)A + ^hi(l + cvA,A) ' ^ 
where f3\ is the second /3-function coefficient. Expanding the right side of Eq. (5.7) up to 
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0(a 3 ) yields: 



a Q = a, [l - a^ \ + a^ 2 A 2 - afaX - a^X 3 + ^aJ/^A 2 + O(aJ)] . (5.8) 



It can be seen that Pi, given by two particle irreducible vacuum polarisation or vertex 
diagrams, occurs only in sub-leading logarithmic terms of the form /3ia™A n_1 . No possible 
re-arrangement of these terms can compensate the manifest gauge dependence of the leading- 
logarithmic terms of the form ((3 a^X) n . 

The property exhibited above, for QCD, of gauge dependence of amplitudes on resum- 
ming one- loop corrections that, at lowest order, are gauge invariant, is expected to be a 
general property of non-abelian gauge theories. In such theories the gauge boson propaga- 
tor, in an arbitary covariant gauge is written as [53]: 



where M is the renormalised gauge boson mass. The topographical structure of diagrams 
contributing to, say, neutrino-neutrino scattering via Z exchange is the same as that for 
quark-quark scattering shown in Fig. 1. The one- loop vertex correction containing the non- 
abelian ZW + W~ coupling is gauge dependent [54]. Since the gauge dependence cancels at 
lowest order (without resummation) then, just as for QCD, it cannot cancel at any higher 
order in the resummed one-loop amplitude. Indeed a similar conclusion as to the necessity of 
the £ = — 3 gauge in order to obtain an effective charge that satisfies a RGE, reached in this 
paper for QCD, has previously been obtained, for the case of the Weinberg-Salam model, by 
Baulieu and Coqueraux [55]. These authors pointed out that, with the special gauge choice 
(in the notation of the present paper) £ = —3, the renormalisation constant of the Z — 7 
mixing term vanishes, so that, in this case, effective charges satisfying separate (decoupled) 
RGE's may be associated the one-loop resummed photon and Z-boson propagators. It is 
also interesting to note that £ = — 3 is the unique choice of covariant gauge for which the 
photon mass counterterms in the renormalised Lagrangian vanish. The case of W exchange 
was not considered, but (as may be seen by inspection of the relevant formulae given, in 
an arbitary covariant gauge, in Ref. [56]) the choice £ = — 3 results in the vanishing of the 
renormalisation constants associated with the one-loop vertex corrections to both the Z and 
W exchange fermion-fermion scattering amplitudes. As for the QCD case considered in the 
present paper, it is then expected that, only for this special choice of gauge, an effective 
charge satisfying a RGE may be associated with the resummed W propagator. 

The pinch technique, and similar methods to formally shift gauge dependent pieces be- 
tween diagrams, have also been applied to electroweak amplitudes [40,46-49]. Although 
gauge invariant boson proper self energy functions may be defined at one-loop level, any 
resummed higher order amplitude is demonstrably gauge dependent, by the same argument 
as that given above for QCD. So, although for a particular choice of gauge parameter (such 
that the sum of all vertex corrections vanish) resummed W and Z running propagators may 
be defined that satisfy a RGE, they cannot, contrary to the claim of Refs. [40,46-49], be 
so defined in a gauge invariant manner. The diagrammatic inconsistency of these proce- 
dures is made manifest by the inclusion in the modified vector boson self energy function of 
box diagram contributions. If the effective charge is expanded as a perturbation series, the 
box diagram contributions at each order will give terms of a geometric series. There is no 




(5.9) 



q 2 -M 2 



,2 _ £ M 2 
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way that such a series can be meaningfully interpreted in terms of a sum of such diagrams 
required by quantum mechanical superposition. In fact, the contributions to physical ampli- 
tudes of box diagrams can be systematically resummed [60, 61, 62], but the result found is 
typically the exponential of a double logarithm of the relevant physical scale, not the sum of 
a geometric series. For the case of the fermion-fermion scattering amplitude the contribution 
of box diagrams is expected to be important only in the \t\ — > limit and to vanish [56] in 
the \t\ — > oo limit. 

A discussion of the gauge dependence, beyond one-loop order, of the resonant Z boson 
amplitude, may be found in Ref. [63]. 

The limited convergence domain imposed by requiring finiteness of the geometric se- 
ries, which occurs in all theories in which the RCC decreases with increasing scales, can be 
avoided by choosing a very high renormalisation scale 8 . This is equivalent, for such theories, 
to the choice, in QED, of on-shell renormalisation, yielding a RCC that is convergent for 
all scales below the Landau scale [20]. Although such a choice guarantees convergence for 
all physical scales below the the chosen renormalisation point, it appears artificial from a 
physical viewpoint. If a strong interaction process at, say, the scale of the mass of the charm 
quark is to be described using a renormalisation point at the GUT scale Qgvt, the formula 
for the RCC at scale m c will depend upon the masses of all strongly interacting elementary 
particles below Qgut and above m c . There will be a phenomenon of 'inverse decoupling' 
whereby the lower the scale the more high mass particles must be taken into account. Feyn- 
man amplitudes using such a renormalisation scale would loose their corresponence (valid 
in the on-shell scheme) with space-time processes. With on-shell renormalisation, the de- 
coupling of heavy particles at low scales is understood in terms of a natural hierachy of 
physical scales. It seems reasonable that the physics of the strong interaction at the scale of 
the charm quark mass, m c , should be independent of the value of the top quark mass, m t , 
when m t ^> m c . According to the the Uncertainty Principle, the contribution to vacuum 
polarisation loops of particles with masses much greater than the propagator virtuality are 
expected to correspond to short lifetime fluctuations giving only a small contribution to the 
radiative correction. This is no longer the case if the RCC is renormalised at scales 3> m t . 

An enormously successful phenomenology of the strong interaction has been developed 
over the past three decades based on perturbative QCD. Many aspects of QCD, such as: the 
existence of the colour quantum number, of spin one gluons with self-coupling, the predicted 
values of colour factors, and a coupling constant that decreases with increasing scales in the 
experimentally accessible region, are confirmed, beyond doubt, experimentally [64]. However, 
it might be hoped that the physical predictions of a candidate gauge theory would be gauge 
invariant at all orders in perturbation theory, as is the case in QED. Explicit calculation for 
the current non-abelian gauge theories (both QCD and electro- weak theory), as reviewed in 
the present paper, seems to show, however, that this is not the case. The point with error bars 
at Q = 90 GeV on the loop gauge curve in Fig. 3 shows the uncertainty on a s (±0.005) of an 
early measurement using hadronic Z decays at LEP [65]. The input value a s (5 GeV) = 0.2 
has been chosen to be consistent with deep inelastic scattering measurements [66]. If the 
asymptotic running coupling constant with time-like argument measured at LEP at the scale 
Q = Mz has a similar value to ofj* \Q) for space-like argument considered in this paper, it 
is clear from Fig. 3 that the measured value of the gauge parameter must be close to —3. 

8 I am indebted to W. Beenakker for this remark. 
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On a more positive note, the phenomeno logical success, in QCD, of calculations based on 
the renormalisation group, as well as of 'renormalon' based models, both of which are shown 
here to require the use of the £ = — 3 gauge, suggests that further understanding requires 
a deeper theoretical explanation of nature's apparent choice of this gauge (see Fig. 3). The 
question is, why, in non-abelian theories, do UV divergent loop diagrams apparently acquire 
logarithmic corrections after renormalisation but not similar vertex diagrams? Since the non- 
abelian triple gluon vertex occurs in both loop (vacuum polarisation) and vertex insertions 
it seems that consistency with experiment requires the cancellation of vertex contributions 
when summed over the complete particle content of the theory. That is, in some conjectured 
modified version of QCD, complete cancellation of triangle vertex insertions similar to the 
anomaly cancellation provided by the particle content of the fermion families of the Stan- 
dard Model. However, consistency with the diagrammatic description (i.e. with quantum 
mechanical superposition) must always limit the scale-range of applicability of the RCC due 
to the convergence properties of geometric series. 

Finally, two important caveats concerning the work presented in this paper should be 
mentioned. Firstly, only covariant gauges specified by the parameter £ are considered, 
whereas certain non-covariant gauges [67] such as axial gauge (A3 = 0) or light-cone gauge 
(n^A^ = 0, n 2 = 0) are freqently employed in QCD phenomenology. The second concerns the 
recently developed successful application of Analytical Perturbation Theory (APT) to QCD 
phenomenology [68] . In order to avoid infra-red divergences when A — > in the conventional 
QCD perturbation series for the RCC, APT introduces, by way of the Kallen-Lehmann spec- 
tral representation of the dressed gluon propagator, the condition of analyticity in the Q 2 
variable. This also imposes a causality requirement. In this case, the correspondence be- 
tween the so-obtained 'Euclidean running couplant', and the sum of a QCD perturbation 
series, in which the terms represent the contribution of specific Feynman diagrams, breaks 
down. No statements may then be made concerning the gauge dependence and convergence 
properties of «e. 
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7 Appendix 

The sum of the first n terms of a geometric series with a negative common ratio r — — \r\ is 
given by the relation [24] : 

S n = 1 - |r| + \r\ 2 - ... + (-[rl)^ 1 = * ~ ^P" - (A.l) 
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It follows that: 



1 + r 



n 



(ra-odd), (A.2) 



S n = ^ -j-q ' (ra — even). (A. 3) 

If |r| < 1 , than = (limit as n — > oo of = 1/(1 + |r|). If |r| = 1 then S* n = 1 for n 
odd and S n — for n even. If |r| > 1, = +oo for n odd and = — oo for n even. 

Table 4 presents values of 5 n versus n, demonstrating the convergence of the series for 
|r| = 1/2 and its divergence for |r| = 2. The Dyson sum of vacuum polarisation insertions 
in the gluon propagator in QCD gives a geometric series similar to (A.l) above. The impos- 
siblity of 'asymptotic freedom', which conjectures that S n — > as n — > oo is evident from 
inspection of (A.2) and (A. 3) above and Table 4. 
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Tables 



i 


Diag 


ram in Fig. 1 


Type 


(ct%/n) 1 x Correction factor, A i: to M w 


1 


b) 


+ 13 <-> 24 


abelian vertex 


%ln^ 

2 u 


2 


c) 


+ 13 ^ 24 


non-abelian vertex 


2 fj, 

^ln« 


3 




d) 


quark loops 


4 


e) 


+ f) + g) 


gluon and ghost loops 


11 



Table 1: UV divergent (before renormalisation) one- loop virtual corrections, at leading log- 
arithmic accuracy, to the 90° quark-quark scattering amplitude in Feynman gauge. 



Gauge 


Feynman (£ = 1) 


Landau (£ = 0) 


Loop (f = -3) 


Vertex (f = 19/9) 


Q L (GeV) 


7.68xl0 8 


1.02x10 s 


301 


oo 



Table 2: Values of the Landau scale Ql (convergence limit) of the QCD effective charge for 
fi = 5 GeV, < = 0.2, n f = 5 



Gauge 


Feynman (£ = 1) 


Landau (£ = 0) 


Loop (f = -3) 


Vertex (f = 19/9) 


Qo(TeV) 


176 


Undefined 


Undefined 


18.1 



Table 3: Values of the Q (the scale at which ot e Jf(Q) vanishes ). Parameters as in Table 2. 
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7(i) (2) (iVLO) 


7(2) (2) (NNLO) 


Gauge 






CfC^ C F C A N f C F Nf 


$ = 1 
Ref. [42] 


13.9 


64 
27 


117.70 -38.50 -|f| 


f = "3 
Ref. [38] 


11.3 




-76.0 23.0 


e = o 

Ref. [38] 


7.6 




-13.2 12.4 



Table 4: Contributions to the n = 2 moment of the non-singlet anomalous dimension in 
deep inelastic scattering on a nucleon due to different classes of Feynman diagrams and for 
different choices of covariant gauge, £. See text for discussion. 



n 


2 


3 


4 


5 


6 


7 


9 


1/(1 + H) 


\r\ = 1/2 


0.5 


0.75 


0.625 


0.688 


0.656 


0.673 


0.664 


0.667 


|r| = 2 


-1 


3 


-5 


11 


-21 


43 


-85 


0.333 



Table 5: Values of the sum, S n , of first n terms of the geometric series in (Al) for |r| = 1/2 
(convergent series) and \r\ — 2 (divergent series). 

Figure Captions 

Fig. 1 Diagrams contributing to Ai^ LO \ Solid lines denote quarks, wavy lines gluons and 
the closed loop in g) ghosts. 

Fig. 2 The topographical structure of diagrams contributing to the resummed quark- 
quark scattering amplitude: a) O(a^), b) O(a^), c) 0(0^). In diagrams containing only one 
vertex insertion the contribution given by the exchange 13 <-> 24 (see Fig. 1) is understood 
to be included. 

Fig. 3 The variation of the Effective Charge (3.5) with the scale Q for different choices 
of the gauge parameter f . 5 GeV < Q < 300 GeV (a e J f (5 GeV) = 0.2, n f = 5). The error 
bars (± 0.005) on the point at Q = 90 GeV on the loop gauge curve are typical of those on 
an a s measurement using hadronic Z decays [65]. 

Fig. 4 The topographical structure of diagrams contributing to: a) the O(a^) con- 
tribution to the anomalous magnetic moment of a quark, b) the four-loop photon proper self 
energy function. 

Fig. 5 The topographical structure of diagrams contributing to non-singlet anomalous 
dimensions in deep inelastic scattering: a) leading order, b) next-to-leading order, c) and d) 
next-to-next-to-leading order. See text for discussion. 
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b) c) d) 




h) i) 

Figure 1: Diagrams contributing to Ai^ ^. Solid lines denote quarks, wavy lines gluons and 
the closed loop in g) ghosts. 
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a) 



b) 




V 2 2 V L L 2 



c) 




V 2 L 2 V L 2 L 3 



Figure 2: The topographical structure of diagrams contributing to the resummed quark- 
quark scattering amplitude: a) 0(a 2 s ), b) 0(a 3 s ), c) 0(aj). In diagrams containing only one 
vertex insertion the contribution given by the exchange 13 <-> 24 (see Fig. 1) is understood 
to be included. 
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0.2 
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Loop (£= -3) 




Landau (£= 0) 




Feynman (c= 1) 




■ ■ ■ Vertex (£= 19/9) 
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i 
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Q (GeV) 

Figure 3: The variation of the Effective Charge (3.5) with the scale Q for different choices 
of the gauge parameter f . 5 GeV < Q < 300 GeV (a e J / (5 GeV) = 0.2, n f = 5). The error 
bars (± 0.005) on the point at Q = 90 GeV on the loop gauge curve are typical of those on 
an a s measurement using hadronic Z decays [65]. 
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a) 




b) 

Figure 4: The topographical structure of diagrams contributing to: a) the 0(af) contribution 
to the anomalous magnetic moment of a quark, b) the four-loop photon proper self energy 
function. 
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a) 0(oc s ) 




b) 



C 2 



c c 



c c 

F A 



0(a 2 ) L+2V gauge invariant 



c) 



C 3 



C 2 C 

^F^A 



C 2 C 

^F^A 



C 2 C C 2 N 

^F^A' ^F^f 



0(a 3 ) L+2V gauge invariant 



d) 



C C 2 

^F^A 



c r 2 c c n 

^F^A' ^F^A^f 



CpC^, CpC^Nj CpC^, CpC^Np CpNj 



0(a 3 ) V 2 +2VL+L 2 gauge dependent 



Figure 5: The topographical structure of diagrams contributing to non-singlet anomalous 
dimensions in deep inelastic scattering: a) leading order, b) next-to-leading order, c) and d) 
next-to-next-to-leading order. See text for discussion. 
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